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ABSTRACT 

Planet-planet collisions are expected during the early stages of the formation of extra-solar planets, 
and are also possible in mature planetary systems through secular planet-planet perturbations. We 
investigate the electromagnetic signals accompanied with these planetary collisions and their event rate, 
and explore the possibility of directly detecting such events. A typical Earth- Jupiter collision would give 
rise to a prompt EUV-soft-X-ray flash lasting for hours and a bright IR afterglow lasting for thousands 
of years. It is suggested that with the current and forthcoming observational technology and facilities, 
some of these collisional flashes or the post-collision remnants could be discovered. 

Subject headings: Planetary systems 



1. INTRODUCTION 

More than 100 extra-solar planets have been detected. 
The growing sample allows us to study the properties 
of these extra-solar planetary systems, and sheds light 
on the origin of our own solar system. At the same 
time, our understanding of astrophysical phenomena has 
been greatly boosted through studying cataclysmic tran- 
sient events such as supernovae. X-ray bursts, gamma-ray 
bursts and their afterglows. Here we discuss another type 
of electromagnetic transient events that arises from colli- 
sions of extra-solar planets. Such giant impacts have been 
long predicted within the theories of planet formation (e.g. 
Wetherill 1990). Giant impacts have been invoked to in- 
terpret the formation of the Earth-Moon (Wetherill 1985; 
Benz, Slattery & Cameron 1986; Canup & Asphaug 2001) 
and the Pluto-Charon (Stern 1992 and references therein) 
systems, as well as the anomalous obliquities of Uranus 
and some other solar giant planets (e.g. Slattery, Benz 
& Cameron 1992). Collisions are also expected in ma- 
ture planetary systems due to slow onset, secular dynami- 
cal instabilities, although these collisions are less frequent 
than in dynamically young systems. The probability of 
collisions is enhanced in white dwarf planetary systems 
because of the perturbation caused by drastic mass loss 
during the formation of the white dwarfs (Debes & Sig- 
urdsson 2002). Accompanying an impact event, a prompt 
transient electromagnetic signal and a fading afterglow are 
expected, but these signals have so far been sparsely stud- 
ied. Even before the discovery of planetary systems around 
main sequence stars, Stern (1994) discussed the frequency 
and infrared (IR) signals of planetary collisions. In this 
Letter, we will sketch the likely electromagnetic signals 
from planetary collisions (§2). We will estimate the possi- 
bility of directly detecting such collisional events and their 
remnants, and discuss the observational strategies (§3). 

2. EUV-SOFT-X-RAY FLASHES AND IR AFTERGLOWS 

As an example, we consider a Jovian planet (primary) 
with mass Mi = 10^° g mj and radius i?i = 7.0 x 
10^ cm rj being hit by a terrestrial planet (secondary) with 

^We note that for collisions between planets with similar masses, 
rather than the signals discussed in this Letter. 



mass M2 = 6 X 10 ' g and radius R2 = 6.4x10° cm , 
where mj and rj are the mass and radius in units of the 
Jupiter values, while and are the mass and radius 
in units of the Earth values, respectively. The reasons to 
choose a Jupiter-sized object being the primary are several 
folds. (1) Over 100 extrasolar giant planets (EGPs) have 
been detected, so that we know they are definitely there; 
(2) The signal is the strongest for these impacts because 
of their large masses (and hence higher Eddington lumi- 
nosities); (3) Because of their larger gravitational potential 
and larger cross-sectional area, the chance of collision per 
planet is higher at higher masses. Since the main emission 
quantities scale with masses and radii, the discussion below 
can be straightforwardly generalized to collisions between 
planets with other parameters (e.g. EGPs). 

The probability of a glacing event is ~ R2/R1 for a 
gravitational-focused system, or ~ {R2/Rif when the 
colliding velocity largely exceeds the escaping velovity of 
the primary. For our parameters, the case is marginal 
for both regimes, and the probability for glacing is be- 
tween ~ 9%(r0/rj) and ^ l%(r0/rj)^, which is neg- 
ligible. Hereafter we will assume head-on collisions 
throughout^. Assuming zero velocity at infinity, the to- 
tal energy of the collision is E' ~ GM2M1/ Ri ~ 6 x 
IQ-io gj-g (^ui^uij/rj), and the relative velocity upon im- 
pact is w ~ 45 km s~^ {'mj/'fjy^'^. Earthquake data re- 
veal that the average p-wave speed interior of the earth is 
~ 10 km (Dziewonski & Anderson 1981), so we can 
generally define a sound speed Cg ~ 10 km s^^Cs, with 
Cs 1 for terrestrial planets. Since v ^ Cs, strong shocks 
will propagate into both colliders, converting most of the 
kinetic energy into internal energy and accelerating non- 
thermal electrons. The shock crossing time for the im- 
pactor, i.e., the terrestrial planet in our case, is 

Tsh ~ R2ICS - 640 s (re/c,). (1) 

This defines the rising timescale of the collision-induced 
electromagnetic fiash. The energy deposit rate during this 
stage \s E ^ E/tsh 9 x 10^"^ erg s~^(m®TO,//rj)(r®/cs). 

chance of glacing events is much larger, which leads to warm planets 



the 
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This is ^ 1.8 X 10'' times the Eddington luminosity 
of the Jovian planet, L^dd = AncGrnpAI/ax ^ 5 x 
lO^'^erg s~^mj (where c, G, nip, ot are the speed of light, 
gravitational constant, proton mass, and Thomson cross 
section, respectively). About one half of the coUisional en- 
ergy is converted to internal energy, of which only a small 
fraction is radiated promptly, while a large fraction is sunk 
as latent heat. The other half of the energy is initially in 
kinetic form, driving a radiation-supported expanding en- 
velope and inducing oscillation and convection within the 
giant planet interior. This kinetic energy will be thermal- 
ized later and radiated away over a longer period of time 
along with the latent heat, as an IR afterglow (see below). 
In the prompt phase, photons must be trapped by the high 
opacity in the regions where the heat is initially deposited, 
and the peak bolometric luminosity is 

Z/pfe = ryZ/Edd ~ 5 X 10^^ erg s"^ r]mj ~ \2>LQr\mj, (2) 

where Lq = 3.8 x lO^'^ erg s^^ is the solar luminosity, ry < 1 
is a factor correcting for less than optimal emission geom- 
etry and other radiative inefficiencies. Since E ';$> L^dd, 
77 ~ 1 is likely achieved. The spectrum is essentially ther- 
mal from a hot spot^ with radius ^ R2 , with a non-thermal 
hard tail causc;d by Comptonization of the thermal pho- 
tons by the non-thermal electrons accelerated from the 
shocks. The typical thermal temperature at the peak time 
is 

= (4^) ^ ^ ^ v'/'my%'/\ (3) 

which is insensitive to both the planet mass and the un- 
known radiation efficiency (i.e. oc rj^^^m^/^). The peak 
frequency is Vpk — 6.6 x 10^^ Hz ry^/^'my^rgj^''^ (wave- 
length ^ 450 A and energy ^ 28 eV). This is in the extreme 
ultraviolet (EUV) band. The peak flux is 

(4) 

where h, k, and a are the Planck's, Boltzmann's and 
Stefan-Boltzmann constants, respectively. The hard non- 
thermal tail is expected to extend into the soft-X-ray band. 

For the thermal spectrum, the flux around v could be 
estimated as F{i^) - {R/DfB^{T)v, where B^{T) is 
the Plank function, R is the emission radius, and D is 
the distance to Earth. One can define a ratio f{v) = 
F{i',p)/ F{v,*) as the fiux contrast between the planet 
and the star at a certain frequency where F{iy,p) and 
F{i', *) are the fluxes of the planet and the star in the 
band around v, respectively. At the flash peak time, the 
EUV flare can greatly out-shine a main-sequence host star 
later than B5, or a white dwarf host star. In the opti- 
cal and IR bands, high contrasts arc also achievable for 
late-type host stars. For example, for a G2 (solar-type) 
star, the U, V, and I bands contrasts at the peak time 
are /(U,pk) - 0.2, /(V,pk) - 0.02, and /(I,pk) - 0.008, 
respectively. Even higher contrasts are achievable when a 
cooler star is adopted. For / > 0.01, the flare is detectable 

^The hot spot will eventually expand via a subsonic excavation flow ( 
is ~ 1 day, much longer than the flash timescale discussed here (eq. [5] ) . 



with photometric monitoring of ground-based telescopes 
(e.g. Borucki ct al. 2001). 

Fixing a particular band, a lightcurve could be depicted. 
The rising lightcurve is very steep, with a time scale of 
~ Tsh ~ 10 minutes. The decay time scale after the peak 
is likely not very long, based on existing observational data 
of other related phenomena, e.g. Comet Shoemaker-Levy 
9 vs. Jupiter collision (Orton et al. 1995), X-ray bursts 
in accreting neutron stars (Lewin, van Paradijs & Taam 
1993) and soft gamma-ray giant flares from magnetars (Fc- 
roci et al. 2001). The prompt radiation pressure would 
drive a fraction of the impacting materials outwards, so 
that initially the Eddington-limited photosphere expands 
and then cools. As the radiation pressure gradually drops, 
the photosphere contracts and the materials falls back onto 
the Jovian planet surface. This fallback refreshes the sur- 
face heating, leading to continually refreshed electromag- 
netic flashes with descending magnitudes. A best guess is 
that the prompt flash will follow a decaying tail extend- 
ing up to several lO's of the rising time scale. The typical 
timescale for the prompt flash during which the emission 
luminosity is above, say 50% of Lpk could be estimated 



n ~ lOTsh ~ 2 hr (re/cs), (5) 



with the averaged bolometric luminosity L^^^^ ~ 
lOLQrjmj. This is the first important time scale in the 

problem. 

The total radiated energy during this prompt flash stage 
is £;/'a'5'i _ Z(_J"'''ri - 3 X lO'*^ erg, which is much less 
than the total impact energy, i.e. <C -E ~ 6 x 10^*^ ergs. 
So the bulk of the collision-deposited energy is not ra- 
diated promptly. Rather, it is mixed into deep inside 
of the planet and stored as heat. It is released during 
a much longer time scale through a decaying thermal- 
radiation afterglow. Model simulations reveal that at later 
times a Jovian planet follows an empirical bolometric lu- 
minosity law (Black 1980; Saumon et al. 1996), i.e., 
L/Lq ~ 10-5m2-35t-i-22^ and Tcf ~ 860 K mY^t^'^-'^'^ , 
where tg is the age in unit of Myr, and the relations are 
valid for > 1. Notice that the luminosity decaying in- 
dex ~ 1.22 > 1, so that the total energy released during 
t > 10^ yr is negligible compared with the energy released 
during the early epoch. At earlier times, generally we ex- 
pect a luminosity decaying law L = Lo{t/to)~°' with a < 1, 
where and Lq are the time and luminosity when such 
a steady fading phase begins. The duration of the bright 
afterglow (the second timescale in the problem) can be es- 
timated as T2 = Tag = [(1 - a)E/Loto]^^'^^-''Ho. The val- 
ues of to, Lq and a, are unknown, and a may also change 
during the decay. In principle, these should be calculated 
from the giant planet models. As a crude estimate, we 
take (tQ.Lo) = (1 yr.3.6 x 10^^L.;.j). The Lq value corre- 
sponds to an initial temperature of Tq ~2500 K (e.g. Stern 
1994 and references therein) emitting from the full giant 
planet surface with the normalized radius rj = 1. Given 
such an initial condition, a < 1/4 is required in order to 
have L{tQ = 1) > 10~^Lq. The bright afterglow duration 

elosh 1989) or a convective flow, but the timescale for this to happen 
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is therefore 



1.4 X 10^ yr {mQmj)r/ 



Tn 



2500 K 



<r2<10Vr, (6) 



where the lower and the higher hmits correspond to o = 

and a = 1/4, respectively, and for the lower limit, the 
explicit dependences on the model parameters have been 
introduced. The peak of this bright afterglow is in the 
IR band, and it is more favorably discovered in longer 
wavelengths. For example, for a G2 host star, the typical 
I- and K-band planet-to-star flux contrasts are /(I, ag) ~ 
2.6 X 10-4 and /(K, ag) ~ 1.7 x IQ-^ (see also Stern 1994). 

Between the prompt flash and the onset of the later 
shallow-decay long-term afterglow, the bolometric himi- 
nosity has to drop by 5 orders of magnitude with a decay 
index a > 1. This would nicely match the (to,io) initial 
condition adopted in the above discussion. In the above 
estimate we have assumed that the prompt flash radiation 
is released from the hot spot, while the long-term afterglow 
emission is emitted from the full surface. It is inferred that 
the Jupiter interior must be convective (Hubbard 1968; 
Bishop & DeMarcus 1970). Thus the collision-deposited 
energy should be rapidly transported to the planet bulk 
during tens or hundreds of the convective overturn time 
scale. The depth into which the impactor velocity is re- 
duced to 1/2 is ~ 2.5 X 10^ cm, i.e., the depth where the 
impacted mass equal to the mass of the impactor, which is 
about 1/3 of the Jovian planet radius^. This length scale 
may be adopted as the convective mixing-length. The con- 
vective speed is less clear due to the unknown interior vis- 
cosity. Nonetheless it should be subsonic, and we take 
Vc ^ 1 km s~^%]c.a for the following discussions. This gives 
a convective overturn timescale ^conv 

0.3?;-o days. The 

timescale for the excavative flows (Melosh 1989), if any, 
is longer than this. Taking a typical overturn number of 
A^c = 100iVc,2 to achieve a full-surface isothermal state, we 
estimate a third timescale 



~ l^Jy at a distance of 10 kpc. Only very nearby events 
(e.g. within 1 kpc) are favorable for radio observations. 

3. DETECTABILITY AND SEARCHING STRATEGY 

Numerical simulations (Chambers, Wetherill & Boss 
1996; Ford, Havlickova & Rasio 2001) suggest that col- 
lisions of the type we are discussing are plausible. The 
basic picture is that there are secular perturbations of the 
inner planets, which over time scales comparable to the 
age of the system (Quinlan & Tremaine 1992) lead to large 
changes in eccentricity and semi-major axis for one or more 
planets, leading to a large probability of collision. This is 
even favorable for those EGPs that are found much closer 
to the host star than Jupiter. Current theories conjecture 
that these EGPs migrate towards the host star within the 
gas disk from which they were born (Lin, Bodenheimer & 
Richardson 1996). The disk dissappears in about 10 mil- 
lion years. Recent meteorite isotope studies suggest that 
terrestrial planets form in a mean time scale of 10 mil- 
lion years (Yin et al. 2002; Kleine et al. 2002; Jacobsen 
2003), so that they could be available for the collisions we 
conjecture''. 

Considering an ensemble of stars with an average age i, 
in order to detect one event with duration r (e.g. n for 
the prompt flash and T2 for the afterglow) after a contin- 
uous observation time of tabs , the critical number of stars 
in this ensemble that have to be searched is 



= ifpNcY 



t 



max(T, to, 



bs) 



(8) 



T3 1 month NcfiV^fl ■ 



(7) 



Since the typical spin period of giant planets is Pp < 1 

day <C T3, one expects a periodic luminosity modulation 
on the EUV-Soft-X-ray lightcurve, caused by the expanded 
hot spot (about 3% of the full surface) entering and leav- 
ing the line-of-sight, over a time scale of T3. This may be 
also detectable in the U, B bands or even in the V band 
through photometry. 

During the vigorous convective epoch, it is likely that 
the convection-induced turbulence would amplify the mag- 
netic field to higher values than conventional (e.g. sub- 
Gauss for Jupiter). The fast relativistic electrons accel- 
erated from the shocks or from magnetic reconnections 
would radiate coherently, giving rise to a radio flare simi- 
lar to the Type-II solar radio flares (but with much higher 
amplitudes). For a flare with brightness temperature 
Tb > 10^^ K, the flux would out-shine the host star, and is 

^The penetration depth at which most of energy is deposited is even deeper, especially when the mantle of the terrestrial planet is stripped 
(Canup & Asphaug 2001), but not much larger than the value quoted here. 

*In case of a late formation epoch for terrestrial planets (Wetherill 1985), a proto-terrestrial planet gives a smaller total impact energy and 
a shorter IR afterglow, but the prompt, Eddington-limited, flash luminosity essentially remains the same (eq.[2]). 

^Although the lack of a major planet in the asteroid belt and the relatively small mass of Mars may be attributed to Jupiter's perturbations 
that inhibited the growth of a terrestrial planet near Jupiter required for collision, we note that the solaj: system might be atypically stable (as 
compared with other extrasolar systems discovered) for anthropic reasons. 



where fp is the fraction of the stars in the ensemble that 
have planets, and Nc is the average total number of colli- 
sions during the lifetime of a typical star in the ensemble. 
Evidence suggests that fp could be close to unity (e.g. 
Lineweaver & Grether 2003; Bary, Weintraub & Kastner 
2002), and we will assume this most optimistic case. Nc 
could be roughly estimated as the number of planets in 
the dynamical system (which is greater at earlier epochs 
since each dynamical interaction destroys one planet) mul- 
tiplied by the probability of direct collisions (rather than 
dynamical ejections), which is in the 40% to 80% range 
(Ford et al. 2001). The high eccentricity of some EGP 
orbits may be caused by dynamical ejection of another 
EGP. Although no collision is involved in such events, the 
eccentric orbits they result more likely lead to strong sec- 
ular perturbations of the remaining planets on a longer 
time scale, particularly so if we assume inclinations are 
increased as well. These actually increase the probability 
of future collisions with planets that might otherwise be 
unperturbed. A good portion of collisions should happen 
in dynamically young systems (which associate with star 
forming regions), but there could be also some collisions 
occuring in mature planetary systems. For an old plane- 
tary system like our Solar system, a typical value of A^c ~ 5 
may be reasonable^, while for dynamically young systems. 
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one may adopt Nc ~ 3. A large fraction of stars arc in 
binaries or multiples. Many of these are either detacihecl 
or very compact, allowing for the presence of planets or- 
biting one or both stars (Holman & Wiegert 1999; Hatzes 
et al 2003). Planets may be less common around stellar 
binaries, but those systems may contribute disproportion- 
ately to the net collision rate, as a stellar companion can 
drive long term secular instabilties among multiple planet 
systems, increasing the collision rate per system. 

Besides A^*, one can also define a characteristic flux f,^^c 
of the coUisional events. Given a number density of the 
stars with the average age t, one can estimate a critical 
distance one has to search in order to find one collisional 
event, i.e., Dc (A^*/n,)^/^. The typical flux can be then 
estimated with D,.- Although A^, is sensitive to the average 
age t of the ensemble of stars being investigated, Dc and 
fv,c are essentially independent of the ensemble adopted, 
given a constant birth rate of stars. 

For the prompt flashes, since there are considerably 
fewer point sources in the EUV and soft X-ray sky and 
since the collisional events out-shine the host stars, one 
does not need to monitor the individual stars in order to 
identify a collisional event. The characteristic flux fi, c is 
relevant to discuss the detectability. The Extreme Ultravi- 
olet Explorer (EUVE) performed a one-year EUV all-sky 
survey (Bowyer et al. 1994) and subsequently ccmducted 
the serendipitous observations with the Right Angle Pro- 
gram (RAP) (McDonald et al. 1994; Christian et al. 1999; 
Christian 2002) for 8 more years. The RAP field of view is 
~ 20 square degrees. So the combined observation is equiv- 
alent to a one-year monitoring with 9 x (20/720) = 1/4 
full-sky coverage. The total number of stars with t ~ 10^" 
yr, for which one collisional event is expected to have oc- 
curred in the EUVE Held of view, is ~ 8 x 10^. The typical 
distance of such an event is ~ 10 — 20 kpc, so that the crit- 
ical peak flux is f,y^c ~ (15 — 60)/xJy peaking around 450A. 
This is the band for which the ISM opacity is the largest. 
At shorter wavelengths where opacity is much lower, e.g., 
50-lOOA (~ 0.1 - 0.2 keV), the flux does not drop rapidly 
due to the non-thermal tail. Even for a flux two orders 
of magnitude lower than the peak flux (a very conserva- 
tive estimate), i.e., (0.15 — 0.6)/L(Jy, this corresponds to 
(170 - 680) coimts ks'^ at lOOA for EUVE (Bowyer et al. 
1994), which is well above the survey sensitivity threshold 
(McDonald et al. 1994). So at least one flare event with 
such a flux level should have been detected. Since during 
~ 10 yrs, ~ 10 X 10"/(10i75) ~ 500 collisional events 
should have occurred in the Milkyway Galaxy (^ 50 per 
year), about 10 more events with lower flux levels might 
have also been caught by EUVE, although with less sig- 
nificance. The predicted flux level is also well above the 
sensitivity of soft X-ray detectors, such as ROSAT, Chan- 
dra, and XMM-Newton, so that there might be such events 
recorded in their archival data as well. A future dedicated 
wide field detector sensitive to (50-200)A would be able to 
detect lO's of EUV-soft-X-ray flares per year. The great 
opacity in the star forming regions would reduce the flux 
level and even the detectable event rate in the softer bands, 
but some events should be directly detected, especially for 
those collisions outside of the star forming regions. 

Another approach to identify prompt flashes is through 
photometrically monitoring a huge number of stars over a 



long period of time. As discussed above, during the flare, 
a 2-20% percent increase of the fluxes are expected for 
the late type stars in the V band and bluer, which could 
be easily detected. Such a photometric search could be 
achieved by future missions such as GAIA (e.g. Perryman 
et al. 2001). Close monitoring ^ 10*^ stars for several years 
should lead to detection of several collisional events. 

The event rate of the IR afterglows is much larger thanks 
to their much longer durations. Due to small flux contrasts 
(^ lO"'"* in K-band), one has to monitor the stars individu- 
ally. The most economic way is to search in the star form- 
ing regions (Stern 1994). For an ensemble of stars with 
f ~ 5 X 10^ yr, one IR afterglow remnant would be found by 
searching for N^^ ~ 1.7 x W f-^{Ncli)-^{i/^ x 10^)T2;i 
young stars, where r2,4 = r2/10^ yr. The signals could 
be identifled through unresolved spectral measurements 
or direct interferometric imaging (Stern 1994). Existing 
groundbased IR interferometers have 1 mas resolutions (J. 
Ge, personal communication), which corresponds to 0.13 
AU at a distance of 130 pc. the distance of the nearest 
star forming region. This is well adequate to resolve bright 
planets from their host stars. 

4. SUMMARY AND DISCUSSION 

We have outlined the electromagnetic signals that would 
accompany with an event of collision between a Jovian 
planet and a terrestrial planet, and discussed the strat- 
egy of possibly detecting these events. The collision itself 
results in a bright EUV-soft-X-ray flare lasting for hours. 
EUVE and some soft X-ray detectors might have detected 
some of these flares. They may be also identified through 
monitoring a huge number of stars with future missions 
such as GAIA, or through future synoptic all-sky surveys. 
After detecting an EUV-soft-X-ray flare, a search of its 
optical eomiterpart through photometric monitoring (like 
catching afterglows in gamma-ray burst study) is desir- 
able. A planetary-rotation-period modulated fading sig- 
nal would be an important clue. Doppler radial velocity 
measurements and IR monitoring may be performed later 
for the collision candidate to verify the existence of the 
planet (s). Each collisional event would leave a bright rem- 
nant glowing in IR. The duration of the afterglow is thou- 
sands of years. Such afterglows could be directly searched 
in the nearby star forming regions, and the planets could 
be directly imaged. The current planet-searching has been 
focused on the planets on stable orbits (Perryman 2000). 
Searching for planetary collisional events will open a new 
window in the planetary science. It allows us to witness 
the gigantic events that are believed to have happened in 
the early history of our solar system, and to directly test 
the planet-moon- formation theories (see also Stern 1994). 
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